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The quantum dynamics of the flat Friedmann-Lemaitre-Robertson- Walker and 
Bianchi I models defined by loop quantum cosmology have recently been translated 
into a spinfoam-like formalism. The construction is facilitated by the presence of 
a massless scalar field which is used as an internal clock. The implicit integration 
over the matter variable leads to a nonlocal spinfoam amplitude. In this paper we 
consider a vacuum Bianchi I universe and show that by choosing an appropriate 
regulator a spinfoam expansion can be obtained without selecting a clock variable 
and that the resulting spinfoam amplitude is local. 

PACS numbers: 04.60.Pp; 04.60.Kz, 98. 



I. INTRODUCTION 

The spinfoam formalism is a covariant approach to quantum gravity closely related 
to canonical loop quantum gravity (LQG) jH|-[9[. There are numerous results on the relation 
between the two languages, for example in the context of the 3d theory [loj], of the full 



hamiltonian theory |11H14|. and in the context of cosmology |15Hl7j. but more clarity is still 
needed. 

The relation between the canonical and the spinfoam languages can be studied in the 
simplified context of cosmology. The canonical quantum dynamics of cosmology is well under 



control in the loop quantization [18H2l|. and a spinfoam expansion has been derived from 



the canonical theory for the flat Friedmann-Lemaitre-Robertson- Walker (FLRW) |22|, |23 
and Bianchi type I [24j cosmologies. In both massless scalar field was included 

in the model. This matter field plays two related roles. First, it allows the theory to be 
straightforwardly deparametrized by treating it as an internal clock variable. Second, it acts 
as a regulator since the implicit integration over it turns distributional transition amplitudes 
into regular functions. This procedure, however, leads to spinfoam amplitudes that are 
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nonlocal in time. A nonlocal spinfoam expansion can still be an effective computational 
tool, but does not match the structure of the spinfoam expression of the general theory 
@, S, |25Hiz|, where locality is a foundational principle and full covariance under the choice 



of clock-time variables is strictly implemented. 

In this paper we study how to obtain a path integral formulation for the loop quantum 
cosmology (LQC) of the vacuum Bianchi I space-time maintaining full covariance under 
the choice of clock-time variables, namely without deparametrizing the theory. We work 
without a matter field and introduce a regulator 5 to control the distributional aspect of 
the transition amplitudes. The physical inner product we define is accurate up to some 



small error which vanishes as 5 — > 0. The idea of such a regulator was introduced in [28 
for the case where the spectrum of the eigenvalues of the Hamiltonian constraint operator 
is discrete, but the regulator used there is not appropriate in the continuous case. In this 
paper we consider two different regulators which we call the Gaussian and the Feynman 
regulators. The choice of the regulator is not trivial: we show that only the second leads to 
local amplitudes. 

The paper is organized as follows: in Section II we briefly review the canonical theory for 



the LQC of vacuum Bianchi I models [2lJ. In Section III we derive two different regulated 



spinfoam expansions for that model and we close with a discussion of our results. 



II. THE LOOP QUANTUM COSMOLOGY OF THE VACUUM 

BIANCHI TYPE I MODEL 

We consider vacuum Bianchi I space-times, possibly with a nonzero cosmological constant 
A. The metric for a Bianchi I space-time is 

ds 2 = -iV 2 dt 2 + ai (t) 2 dxi + a 2 {t) 2 dx 2 2 + a 3 {t) 2 dx 2 3 , (2.1) 

where N is the lapse and a« are the three directional scale factors. In previous treatments 
of the Bianchi I model in LQC, the lapse has been chosen to be N = |aia2fl3| and we make 
the same choice here. 

We introduce a fiducial cell in order to deal with trivial divergences in the canonical 
analysis (the resulting physics is independent of the choice of the cell). We use a cell 
adapted to the symmetries of the space-time: a rectangular cell with lengths £x,£% an d l 3 
with respect to the fiducial metric ds 2 = dx 2 + dx 2 + dxg. 

The variables used in LQC are the Ashtekar connection A 1 and the densitized triad Ef; 
in Bianchi I space-times they can be parametrized as follows [2l| : 

4=£(d*<). and E t ^{£j, (2.2) 

where there is no sum over the i and V Q = l\li£%. The Poisson brackets are given by 

{ci,Pj} = SKGjSij, (2.3) 
and the (c,p) variables are related to the ones in the metric by 

Pi = sgn(ai) I a 2 a 3 \£ 2 £ 3 and c x = * dQl , (2.4) 
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where 7 is the Barbero-Immirzi parameter. The other relations can be obtained via permu- 
tations. 

Following the ideas of LQG, the fundamental variables of the kinematical Hilbert space 
are holonomies of the connection and area fluxes of the densitized triads rather than the 
connection and densitized triads themselves. Accordingly, only operators corresponding to 
complex exponentials of c, are defined in the quantum theory, no q operator exists. Since 
the j>i correspond to the physical areas of the fiducial cell, they can be directly promoted to 
be operators. See [2l[ for the details of the construction of the kinematical Hilbert space. In 
the pj-representation the resulting Hilbert space is composed of functions ip(pi,P2,P3) with 
finite inner product 

H 2 = ^(P1,P2,P3)^(P1,P2,P3)' (2-5) 

On these states Pi act via multiplication and the holonomies e lfJ-iCi act as translations. Finally, 
the kinematical inner product of basis states is 

(PI,P2,P3\p'i,P2,P3) = Wi'WWs* ( 2 ' 6 ) 

note that the <5's are Kronecker delta functions, not Dirac delta distributions. 

The Hamiltonian constraint operator has been derived in 21] . Before writing the explicit 
form of the operator itself, we make a change of coordinates which simplifies its form. We 
introduce the A,- variables defined as 



Pi 



sgn(A i )(47r 7 v / A4i) 2/3 A-, (2.7) 



where Af 2 ^ = Ay/Sw^i^ is the area gap in LQC, the minimum eigenvalue of the area 
operator. Next, let 

v = 2A!A 2 A 3 , (2.8) 

which is proportional to the physical volume of the fiducial cell. 

Any state 1 can now be described by a ket |Ai, A2, v). The kinematical inner product is 

(Ai, A 2 , v\X[, A' 2 , v') = 5 \ 1 \' i 5 \ 2 \' 2 5 vv i . (2.9) 

As shown in }21| . the Hamiltonian constraint operator for the vacuum Bianchi I model 
(with a minor modification in order to include the cosmological constant A) is 2 



v 



+ 2)y/v(v + A) tf+(Ai, A 2 , v) - v(v + 2) +(Ai, A 2 , v) 



16 

- v{v - 2) % (Ai, A 2 , v) + (v- 2)< s /v|«-4|tfj(Ai, A 2 , v) 

2 A A 111 „.2 



+ 4 7 ^AA^*(A 1 ,A 2 ,u) 



(2.10) 



1 This is not true for states corresponding to singular geometries, but since these states decouple under the 
dynamics their behavior is trivial and not important for this work. 

2 Although in general the variables Xi and v can be negative as well as positive, it is possible to only consider 
the positive case due to the parity properties of the wave function, see 21] for details. This is what is 
done here. 
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where the \l/^ 4 are defined as follows: 

^(Ai, A 2 , u = * — — ■ A 1; A 2 , v ± n + * — — ■ A 1; A 2 , v ± n 

\f±2 f / \f±2 

T A; ±2 , v±n A \ _ /v±2 . . 

+ ^ A i>^77 ■ A 2 ,v±n +* Ai,A 2 ,v±n 

\ v v ±2 J \ v 

+ Ai,— — .A 2 ,v±nJ + *f Ai,^^ -A^vinJ . (2.11) 

The Hamiltonian constraint defines the physical inner product. There are two ways 
that have been used to construct the physical inner product in cosmology. The first is 
to "deparametrize" the theory by choosing one of the variables as a time-variable, and 
expressing the dynamics of the other variables with respect to it. This is commonly done in 
loop quantum cosmology by coupling a massless scalar field to the metric and using it as a 
clock. 

The second possibility, which is the one we use here, is to maintain explicit covariance 
under the choice of the independent variable and define the physical physical inner product 
by "group averaging". Heuristically, given a suitable state in the kinematical Hilbert space 
a physical state is obtained by acting on the states with a delta function of the con- 
straint, |$) P hy = 5{Ch)\4>)- More precisely, the physical inner product can be defined by the 
expression 

1 r°° 

WW = tt da (^ aCH |0>kin. (2-12) 
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Denoting (Ai, A 2 , v) by A, the physical inner product between the physical states generated 
from the kinematic basis states for the vacuum Bianchi I model in LQC is given by 

-. -, 1 f°° -> - -> 

(AiHAo> phy = t- da (\ F \e taCH \X ) kiQ . (2.13) 
2tt < /_ 00 

This expression may contain divergences. In the next section we reformulate it so that it 
can be computed explicitly. 

III. THE VERTEX EXPANSION 

In this section we express the physical inner product via a Feynman path integral con- 
struction with the form of a vertex ex pan sion, much like spinfoam models. In the first part, 
we review the procedure presented in [22 -24]| for the FLRW and Bianchi I models and show 



how, when it is naively followed, the final result fails to be well defined for vacuum space- 
times. In the following subsections, we present two examples of how one can introduce an 



external regulator as proposed in [28| in order to obtain a well-defined result. 



A. The Standard Procedure 

As usual in a Feynman path integral construction, we break up the action of Ch into 
equal parts (i.e., exp(iaCn) = [exp(iaCn /N)] N ) and insert a completeness relation between 



5 



each term: 

- - 1 f°° - - 

(AF|A )p hy = ^y da(X N \e iaC "\X ) 



= ^J da^^... ^(A F |e^|A^_ 1 )...(A 2 |e^|X 1 )(A 1 |e^|Ao), (3.1) 

Ai A2 Ajv-i 

where we have introduced e — a/N. Since this relation holds for any N, it is straightforward 
to take the limit of iV — > 00. However, an important point here is that the N — > 00 limit is 
inside the integral over a. 

Now, in the limit of small e (or large N) one immediately sees that (up to higher order 
terms in e which can safely be neglected) 

(A|e ieC "|A) «e iee xx and (X\e iedH |A') « z'e% x „ (3.2) 

where 0?r, = (A|Cf/|A') are the matrix elements of Ch and in the second relation we have 
assumed that A 7^ A'. 

We now reformulate the inner product in terms of a vertex expansion. This is done by 
rewriting the sum in terms of the number M of transitions in A, that is to say the number 
of times Aj+i 7^ Aj. Denoting the 'time'-step of the i-th transition by N iy one finds that (see 
[22I . 23] for a derivation of this result) 

(A F |A ) phy = — / da lirn^ J2 " " " J2^ e Th X F x M J ■ ■ ■ fa^Xo) 



A/=0 iV, 



Am- 



e 



(e^x^^-^i-i^XoXo)^- 1 , (3.3) 
where we have introduced the shorthand 

N-l N M -1 N 3 -1N 2 -1 

E= E E •■•EE- < 3 - 4 > 

N M =M iV M _i=M-l AT 2 =2iVi=l 

In the iV — )■ 00 limit, the sums become integrals and the vertex expansion for the physical 
inner product is given by 

(A F |Ao)p hy = — / da J2 ■ ■ ■ a f a m _ 1 • • • e XiX A &> ■ • • > ^/-i; a), (3.5) 

M=0 =* 



where 



A(A l5 . . . , A M -i; a) =* A/ / dt M / dt A/ _! . . . / dt^e^vv)-^ 



JO 



( e ie A M -i^M-i )*m-*m-i _ _ _ (e^XiXi ( e i6 Vo )*i . (3.6) 
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The integrals in this expression can be evaluated, the general solution is (see J23J for details) 



Am-i! a) 



n 

k=l 



n k -l)\ \oe 



Hfc-l 



fclHWe.-e,)' (37) 



where Bj are the p distinct values of A taken along the history and rij is the number of 
times the value 0j is repeated. In the simplest case of all the ni being 1, this gives 3 



M 



A(Ai, . . . , Am-i! op = 2J 



(3.8) 



=0 



The last step is to perform the integral over a and express (Air|A ) p hy as a sum over M 
of some function of the intermediate steps Aj. However, this requires pulling the infinite 
sum over M outside of the integral and this causes major difficulties: in the simplest case 
where all A, are different, one obtains a sum of Dirac delta distributions, not a function. In 
the more general cases, the situation is even worse as one obtains derivatives of Dirac delta 
distributions as well. While this provides a formal perturbative solution to the constraint, 
it is distributional term by term. This would give a sum over distributions outside of any 
integral and the vertex expansion presented here would fail to be well defined. This procedure 
is carried out in the Appendix for the sake of completeness. 

If a massless scalar field is added then the distributions appear under an integral over 
the scalar field momentum. Similarly, any degree of freedom with a continuous spectrum 
introduces integrals which give meaning to the distributions which arise. It is the totally 
discrete nature of the kinematic Hilbert space in the vacuum case that gives rise to the 
ill-defined expansion above. As there is no matter in this case, we introduce an external 
regulator 5 in order to obtain functions rather than distributions in the expression of the 
physical inner product. There are two natural regularizations which we present here: a 
Gaussian suppression of the integral at infinity and an exponential suppression a la Feynman. 



B. Gaussian Regulator 

Recall that the group averaging procedure can be thought of as putting a Dirac delta 
distribution projector in the kinematical inner product. A natural way to regulate the vertex 
expansion is to approximate the Dirac delta distribution by a Gaussian with a small spread 
5. The approximated physical inner product given by this regulated group averaging is 

(A f |Ao)g, 5 = da (A F |e^- 52a2 |A ) kin . (3.9) 

While the projector given by group averaging kills all of the parts of the given wave function 
that don't satisfy the constraint, this regulated 'projector' exponentially damps all of the 



3 Even though expression ()3.7j) is well defined in the limit that some of the diagonal matrix elements of 
the constraint are zero, the integral over a in Eq. (|3.5j) is divergent instead of distributional in this limit. 
For this reason, we have added a cosmological constant to the model since the constraint has no diagonal 
terms if A = 0. Alternatively, we could have instead considered a different basis in the kinematical Hilbert 
space or worked with the master constraint. 
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parts of the wave function that do not satisfy the constraint. Clearly, the correct projector 
is recovered in the limit 8 — > 0. 

Following the same procedure outlined above in order to obtain a vertex expansion for 
this regulated inner product, we obtain 

-| /"OO °° 

(X F \X ) G , 5 = Y \ daJ2J2 X,a m _ 1 • • • ©^"^(A!, . . . , X M -u a), (3.10) 

where A(Xi, . . . , Xm-i] ot) is again given by Eq. (13. 7ft . Pulling the sum over M outside of 
the integral and evaluating the integral over a, we have 

(AV|Ao> g , 5 = — E e XrX*-i ■ ■ ■ ^XoMs^i, Xm-i), (3.11) 

note that A(Xi, . . . , Xm-i) represents the "group-averaged" A(Xi, . . . ,Xm-i',oc). For the 
Gaussian regulator presented here, one can see that 

p i / Q \ "fc- 1 P e -e?/45 2 

A GA M,...AM- 1 ;a)=U ¥ ^\ae k ) Y t - n ^e,-e 1 y < 3 ' 12 ' 

Again the 0j label the p distinct values of 6^ taken along the history (Ao, . . . , Am) and rij 
the number of times that value is repeated. Notice that one cannot take the limit of 5 — > 
at this point, the fact that S is nonzero is absolutely necessary in order to pull the sum over 
M outside of the integral. 

The integral here gives Gaussians and derivatives thereof instead of Dirac delta distribu- 
tions and derivatives thereof. The regulated physical inner product is then well defined and 
is a good approximation: 

(A^IAo/phy ~ (A_f|Ao)g,<5- (3.13) 

This expression can be used as a computational tool but since it is nonlocal, it cannot be 
viewed as a building block for the local spinfoam expansion of full quantum gravity. 



C. Feynman Regulator 

Another possible approach is to split the integral given in Eq. ( 12 .131) into integrals over 
positive and negative a and regulate each integral separately, much as is done for the Feyn- 
man propagator, with an exponential suppression as follows: 

(Af|Ao)p, 5 = — / da(A F |e^- fo |A ) kin + — / da (A F |e^ +5a |A ) kin (3.14) 

= (Xf\Xo) + ,s+(Xf\X )-,8i (3.15) 

this approximates (Ai?|Ao) p h y up to some small error which vanishes as 5 — > 0. 

As an aside, it is interesting to note that if we construct a phase space path integral for 
this model as done for FLRW in 29] we find that each half is related to fixing one sign for 
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the lapse or fixing a single direction for the time evolution. It is necessary though to include 
both terms to have a solution to the constraint. 

Now, the integral over negative a is the complex conjugate of that over positive a, so 
for simplicity we first focus on just one half of this expression and see what effect the 
regularization has on the integrals over a. Taking 

_ _ i r°° 

(\ F \\o)+,5 = tt da (AF|e^~ 5a |A ) kin , (3.16) 
27r Jo 

the expansion of (Ai?|e \Xo)km can be carried out as before and then the integral over a 
gives 

1 oo 

<aV|Ao> +> , = - J2 E •••E e AVM- 1 "-- e A 1 Xo^(^"--^M-i), (3.17) 
where the amplitude is given by 

POO 

A +)S (Xi,...,Xm-i)= / d«e- 5a A(A 1 ,...,X A f-i;«), (3.18) 



and A(Xi, . . . , Xm-i] oc) is given by Eq. ( 13. 7(1 . 

Surprisingly, the integral over alpha reduces the very nonlocal expression given in Eq. ( 13.7(1 
with a simple local one being just a product of matrix elements: 

i(-l) M 

A + , 5 (A 1; ■ ■ ■ , = n M : n J ~ . r\ ' ( 3 - 19 ) 

Putting all of the pieces together, the expansion for the positive half is 

A/=0X - \- A A ./A A / / V A A / 



w— 1 



The 5 remains non-zero in the resulting expansion, since in the limit that S goes to zero it 
reduces to a sum of distributions. The physical inner product, regulated by the parameter 
5, is given by the sum of this expansion with its complex conjugate. Equivalently, 

(A F |X ) phy = 2K((A f |Ao) +5 ) +err(5), (3.21) 

where err(<5) is a small error term which vanishes in the limit of 5 — > 0. Notice that ( 13 .20(1 
can be rewritten in the form 

oo 

(X F \X ) + s = ± E E UMXf) HA v (X f ), (3.22) 

where / = 0, M; v = 1, M; the "face" amplitude is A/ = (9 X r + itS) -1 , the "vertex" 
amplitude is = —9^ ^ , and we write Xm ='■ Xp. The expression (13.221) is precisely the 

XvX v _\ 

general expression for (local) spinfoams, if we identify the sum over spinfoam two-complexes 
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with the sum over M, the sum over colorings with the sum over the A/, the spinfoam vertices 
with the transitions, and the spinfoam faces with the sequences of steps without transitions. 

The fact that the Feynman regularization gives an amplitude which is local is a remarkable 



result. In the previous works [22H24J. |28| . the reconstruction of the spinfoam-like expansion 
obtained from LQC lacks locality in the resulting expression; here locality means that the 
amplitude of one history is the product of amplitudes of the elements forming the history. 
Notice that in quantum mechanics locality is always present in sums over histories. For 
instance, in quantum electrodynamics Feynman graphs, the amplitude of a graph is the 
product of the amplitudes of the vertices and of the propagators. In the functional integral 
formulation of quantum mechanics, the amplitude is the exponent of an integral, which is 
to say (morally) the product of exponents: again a product of local terms. At least in the 
Feynman graph integral, and possibly in other contexts, a nonlocal measure might indicate 
that something is wrong. Thus, this type of locality is a desired feature of a sum over 
histories formulation of quantum gravity. 

A naive regularization is not enough to obtain locality, as we have seen in the case of 
the Gaussian regulator where one again obtains a nonlocal amplitude. With the Feynman 
regularization, the amplitude is local. 



IV. DISCUSSION 

We have shown that in loop quantum cosmology it is possible to introduce a regulator in 
order to approximate the vertex expansion of the physical inner product and that one such 
regulator gives a local expression. The strategy to calculate the physical inner product is 
to choose a small, nonzero 8 and to perform one of the regulated expansions presented here 
to some number of transitions m and this approximates the true physical inner product to 
some order determined by the regulator 5 and the number of transitions m. A priori, the 
convergence properties of the series is unknown. 



The idea of introducing an external regulator was proposed in [28| in the context where 
the Hamiltonian constraint operator has a discrete spectrum. In that case, the regulator 
removes actual divergences in the expansion. In the case of a continuous spectrum, the 
terms in the expansion are simply distributional. The regularizations introduced here are 
natural from the point of view of handling distributions. 

The Feynman regularization adds an additional similarity to spin foams. The spin foam 
amplitude for a single triangulation and set of labels, for a single history of spin networks, 
is given by a product of amplitudes associated to each vertex, face, and edge. In contrast, 
for the Gaussian regulator presented here and for the vertex expansion of FLRW with a 
massless scalar field given in (22|, 23[ the amplitude for a single discrete history is not a 
simple product of amplitudes, rather it is a nonlocal expression depending on the properties 
of the entire history. By not integrating away the matter degrees of freedom from the bulk, 
and using the Feynman regulator introduced here, the amplitude is reduced to a simple local 
product. 

The local form of the amplitude makes the relation between the expansion of the LQC 
scalar products and the spinfoam formalism transparent: the sum over the number M of 
transitions is recognized as the analog of the sum over two-complexes, the sum over the A/ as 
the sum over spinfoam colorings, the transitions as the spinfoam vertices and the sequences 
of steps without transitions as the spinfoam faces. 
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In addition, we observe that by using the Feynman regulator in the models studied in 22 



24J ] before integrating over the momentum of the scalar field, one obtains an integrand which 



is well-defined on its own (rather than distributional) and which gives a local expression for 
the amplitude, thus providing a local extension to the results of these previous works. 
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APPENDIX: Distributional Expansion 

If we ignore the problems mentioned in Sec. IIII Al and evaluate the integral over a without 
introducing any sort of regulator, the resulting expansion is 



(AiHAo)phys E ■ ■ ■ E A F A M 



_! •••@A 1 A 74 ( A l'---' A ^-ly 



(1) 



where 



A(Ai, . . . , A M _i) = — da A(Xi, . . . , A M -i; a). 

Z7T 



This integral can be easily evaluated giving 



A(Xi, . . . , Am-i) 



n 



L (n k - 1)! \dG 







rifc-1 



E 



(2) 



(3) 



where, as before, 0j are the p distinct values of @X \ taken along the history and Hi is the 
number of times the value 0j is repeated. One can verify that this distribution is equivalent 
to the following simpler one: 



A(A 1 ,...,A M _ 1 ) = (-1) M+1 E 



i=i 



n, - 1)\ yoe 



5(6,) 



rp 9 nj 



(4) 



At a formal level, this sum of distributions provides a perturbative solution to the constraint. 
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